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Abstract. Graphene epitaxially grown on Ru(OOOl) displays a remarkably ordered 

pattern of hills and valleys in Scanning Tunneling Microscopy (STM) images. The 

extent to which the observed "ripples" are structural or electronic in origin have 

^ ■ been much disputed recently. A combination of ultrahigh resolution STM images 

■"'^ I and Helium Atom diffraction data shows that i) the graphene lattice is rotated with 

respect to the lattice of Ru and ii) the structural corrugation as determined from He 
diffraction is substantially smaller (0.15 A) than predicted (1.5 A) or reported from 
I _. ^ X-Ray Diffraction or Low Energy Electron Diffraction. The electronic corrugation, on 

f-^ • the contrary, is strong enough to invert the contrast between hills and valleys above 

^^ , +2.6 V as new, spatially localized electronic states enter the energy window of the 

STM. The large electronic corrugation results in a nanostructured periodic landscape 
of electron and holes pockets. 



PACS numbers: 73.20.-r, 68.37.Ef, 68.55.-a, 81.05.Uw 
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1. Introduction 

The growth of large, highly perfect epitaxial graphene monolayers is a prerequisite for 
most practical applications of this promising material [H 12]. In addition, it is crucial 
to understand the interaction of graphene with the surfaces of substrates of different 
nature (oxides, semiconductors or metals), as well as with adsorbed molecules, in view 
of the sensitivity of the conduction properties of graphene to them [3]. Nanostructuring 
graphene (in stripes, dots or by periodic potentials), in turn, may reveal new physical 
phenomena and fascinating applications |1] . Most of these topics can be characterized in 
detail in what has become a benchmark system for epitaxial graphene: a self-organized, 
millimeter-large, periodically "rippled" epitaxial monolayer of graphene grown by soft 
Chemical Vapor Deposition under Ultra High Vacuum (UHV) conditions on single 
crystal metal substrates with hexagonal symmetry, such as Ru(OOOl) [HI El [7], Ir(lll) 
[H [9] or Pt(lll) [To]. The superb control that allows the UHV environment facilitates 
a precise characterization of the system down to the atomic scale. 

The difference in lattice parameter between graphene (gr) and the different metal 
substrates causes the appearance of moire patterns with a range of apparent vertical 
corrugations and lateral periodicities with respect to the basic graphene structure 
[3 El [71 m m [10]. In-situ STM imaging of graphene monolayers on Ru(OOOl) reveals 
periodic corrugations with (12xl2)[5l [HI [7] and surface X-ray diffraction (XRD) 
experiments give a (25x25) periodicity based on the distortion of the first Ru layer under 
the graphene [TT]. The carbon atoms on the gr/Ru(0001) are electronically inequivalent 
as reflected in the X-Ray Photoelectron Spectroscopy C Is core levels, which show two 
peaks, i.e. two differently bonded C atoms [T2] . 

The apparent amplitude in STM of the corrugation of the ripples for gr/Ru(0001) 
decreases from 1.1 A to 0.5 A when the tunneling bias goes from -0.8 V to 0.8 V [T3] . 
While STM studies [6] ascribed most of the apparent corrugation to electronic effects in a 
weakly (<0.3 A) structurally corrugated graphene overlayer. Density Functional Theory 
(DFT) calculations [T4l IT5] , assuming that the graphene monolayer is strictly aligned 
with the Ru lattice, have concluded that the ripples are mostly structural in origin, 
predicting a geometrical corrugation of 1.5 A. On the contrary, a very recent ab initio 
calculation that assumes a slightly different registry between graphene and Ru(OOOl), 
gives a corrugation value of 0.24 A [16]. A first fit to surface X-ray Diffraction data 
suggest a corrugation of 1.5 A [TT], using the same technique (XRD) a recent fit [T7] gives 
a corrugation of 0.82 A. Finally in a recent work using Low Energy Electron Diffraction 
(LEED) a corrugation value of 1.5 A [T8] have been proposed. It is worth to mention 
the limited sensitivity of diffraction techniques to the position of the light Carbon 
surface atoms as compared with the Ruthenium ones. Recently, a model calculation 
[T9] has found, contrary to the previously mentioned DFT calculations [111115], that the 
structural corrugation of the moire pattern is comparable to the atomic corrugation of 
the carbon atoms. Considering its paradigmatic nature as a prototype nanostructured 
graphene system, it is rather upsetting that there is still a lack of consensus regarding 
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the structure of gr/Ru(0001). 

In this work we study the geometric and electronic structure of graphene epitaxially 
grown under UHV conditions on Ru(OOOl) by low-temperature STM and Helium Atom 
Scattering (HAS). The high-resolution STM images of large, atomically perfect domains 
allow us to determine that the graphene lattice is rotated by 0.5° with respect to the 
Ru lattice. The moire pattern, which is rotated by 5.0° with respect to the Ru lattice 
in the STM images, acts as a magnifying lens for the small angle of rotation between 
the two atomic lattices [20]. The intensity of the He diffraction peaks, allows us to 
quantify accurately the apparent corrugation of the total charge at the external surface 
(at about the same distance explored by STM), i.e. the closest experimental measure 
of the geometric corrugation of the ripples, which turns out to be of the substantially 
smaller than predicted [T4| [T5] or reported by XRD [HI [17] and LEED [18]. The much 
larger apparent STM corrugation depends strongly on the bias voltage, and can even 
be inverted above -1-2.6 V, revealing that, in addition to the geometric corrugation, a 
much stronger electronic corrugation exists. 

2. Results and discussion 

The STM and HAS experiments were carried out in different UHV chambers with 
base pressures of lxlO~^° mbar, equipped with standard facilities for metal surface 
preparation, ion gun and mass spectrometer, gas exposure. Low Energy Electron 
Diffraction and Auger Electron Spectroscopy. Atomically clean, bulk C depleted, 
crystalline Ru(OOOl) surfaces were prepared by standard sputter/anneal procedures 
followed by oxygen exposure at 1150 K and a final flash to 1400 K in UHV, which 
resulted in large terraces displaying atomic resolution, separated by monoatomic steps. 

In order to minimize the possible influence of defects and domain boundaries in 
the observed registry between the graphene overlayer and Ru(OOOl) that may induce 
distortion on the moire superstructure [5], special care was taken to clean the bulk of 
the Ru crystal from C by repeatedly exposing it to 3xl0~^ mbar of O2 at 1150 K and 
flashing it to 1400 K. Different partial pressures of ethylene were explored during the 
growth. The samples of graphene were prepared on multiple occasions on two different 
Ru(OOOl) single crystals following the procedure that yields almost perfect overlayers 
and larger domains: the Ru crystals, kept at 1150 K in Ultra High vacuum (UHV), were 
exposed to ethylene at pressures of 2x10^^ mbar for 3 minutes (48 L, 1 L=1.33xl0^^ 
mbar s) or 3xl0~® mbar for 10 minutes (24 L). These exposures were enough to saturate 
the surface. The temperature was held at 1150 K for further 2 minutes after removing 
the C2H4 gas from the chamber. 

These optimum conditions were determined from the combined STM and HAS 
analysis of the grown graphene films, which display almost atomically perfect domains 
covering completely the surface. Figure 1(a) shows an STM image of a single domain 
with a lateral size of 2000 A. The periodic bright dots are the ripples of the moire pattern 
of epitaxial gr/Ru(0001). Only one domain is visible in the image as reflected in the 
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sharpness and order of the spots in the 2D Fast Fourier Transform (FFT) of the image 
shown on the inset. Note that the ripples in all the terraces are aligned exactly along 
the same direction. The only defects present in the image are bubbles of Argon (green 
circle) buried in the substrate during the ion bombardment of the Ru(OOOl) crystal 
[211 122] that appear as cloudy regions. They do not perturb significantly the periodicity 
of the moire superstructure. 

In larger images, such as the one shown in Figure 1(b), the presence of dislocations in 
the moire structure (black circles) is a clear indication of the existence of several domains. 
In order to elucidate if the domains are due either to the coalescence of graphene islands 
with a different registry respect to the Ru(OOOl) lattice or to misalignments between 
the graphene lattices on the different domains, we have calculated the FFT of the 
image |23]. The result is shown in the inset in Figure 1(b), there is only one set of 
spots in an hexagonal arrangement corresponding to the moire superstructure. From 
the comparison with the FFT shown in Figure 1(a) it is clear that these spots are 
elongated. The origin of this deformation is the presence of different domains with a 
small misalignment between them. The angular spread of the moire patterns can be 
determined by measuring the angle between the green lines that connect the (0,0) point 
with the boundaries of the elongated (1,1) spot and turns out to be 10°. 

Figure 2 (a) shows a high resolution STM image of the graphene monolayer, with its 
characteristic triangular array of bumps separated (29.3 ± 0.8) A and the simultaneously 
resolved atomic C lattice. The green dotted line indicates the high-symmetry [1120] 
direction of the carbon lattice and the blue dotted line the corresponding high symmetry 
direction of the hexagonal moire pattern. The atomic rows are clearly not aligned with 
the ripples in all the areas free from defects that have been examined. The angle 
between both directions is ^Pgr,moire = 4.5° ± 0.5°. Similar misalignment between the 
carbon lattice and the ripples can be observed in LEED patterns measured with a Low 
Energy Electron Microscope (LEEM) on a single domain island of epitaxial gr/Ru(0001) 
but went unnoticed [TJ [23] . For diffraction patterns measured with conventional LEED 
on continuous graphene films the co-existence of the misaligned domains makes the 
diffraction spots wider but centered on the high symmetry directions (see insets in 
Figure 1(b) and in Figure 4(a)) rendering difficult to identify such a small rotation 
between the atomic lattice and the moire superstructure. The characterization done 
in real space by means of atomically resolved STM images published so far show only 
one or two unit cells of the moire superstructure, making difficult to identify such small 
rotation but in the high quality ones [251 126] it is easy to see the rotation between the 
atomic rows and the ripples. This fact, however, has not been noticed. 

In Figure 2 (b) the 2D Fast Fourier Transform of the topographic image shows 
simultaneously the periodicity of the moire superstructure and the carbon lattice (green 
circles). From the power spectrum it is clear the co-existence of two different periodicities 
rotated respect each other and with a distance ratio of 1/12 approximately. Figure 
2 (c) shows the profile in the spectral density drawn in the FFT image along the 
line connecting the spots corresponding to the graphene atomic periodicity i.e. along 



itaxially Gr 




O 



\. 





Figure 1. (a) 2000x2000 A^ STM image of grapliene on Ru(OOOl). The color scale 
has been adjusted independently on every terrace. The ordered white dots correspond 
to the maxima of the moire superstructure. The green circle encloses a subsurface 
Argon bubble in the Ru(OOOl) substrate. The inset shows the Fourier transform of 
the image showing the existence of a single hexagon and well defined spots, a clear 
indication of long range order on the moire superstructure, (b) 3000x3000 A^ STM 
image of graphene on Ru(OOOl). The STM images have been differentiated along the 
horizontal direction. Black circles mark some of the dislocations present on the moire 
superstructure. The inset shows the 2D Fast Fourier Transform of the image containing 
the first hexagon of the moire superstructure. The spot corresponding to the moire 
periodicity is elongated showing an angular distribution of ±5° marked by green lines. 



the high symmetry direction of the substrate (green circles in panel (b)). The peaks 
corresponding to the atomic periodicity can be seen at both ends of the spectral density 
(marked with red arrows). Due to the small rotation between the carbon atomic 
structure and the moire superstructure, only the first two peaks corresponding to the 
ripples periodicity can be seen. On the contrary a line profile in the spectral density 
rotated by 5° pass through all the peaks corresponding to the periodicity of the moire 
superstructure as shown in Figure 2(d) (blue profile). 

The rotation between the C and the moire lattices for gr/Ru(0001) reflects in an 
amplified fashion the misalignment between the graphene monolayer and the underlying 
Ru lattice [20]. In order to determine directly the epitaxial relationship of the graphene 
lattice with the Ru substrate, the growth of graphene has been stopped just below 
the completion of the monolayer [27] allowing the existence of small patches of the Ru 
substrate not yet covered by graphene. Since it has proven very difficult to obtain 
simultaneously atomic resolution in both the graphene islands and the Ru substrate, 
the surface has been exposed to a small O2 dose. Oxygen does not adsorb at 300 K 



ially Gr 




(b) 











^^^B^^^^^l 


"*-^^ 


-^ 


■■|^r^/.v- ''jjl^^^^^^^^l 


^~~~~@ 














"^^^^H^VY 


4' 



-1,2 -0,8 -0,4 0,0 0,4 0,8 1,2 

Wave vector (A"') 




Wavevector (A" ) 



Figure 2. (a) Atomically resolved STM image (130x 130 A^) of tiie complete graphene 
monolayer grown on Ru(OOOl) recorded at 4.6 K in an area free of defects (image 
setpoint V^ = +lmV and It = 1 nA). The green dotted line in the topographic 
image follows the high-symmetry [1120] direction of the C atomic rows in the graphene 
layer. The blue dotted line shows the direction of the moire superstructure, (b) Power 
spectrum of the STM image shown in (a). The green circles mark the spots due to the 
graphene atomic periodicity; (c) Spectral density along the high symmetry direction 
of the atomic structure marked with a green line on panel (b). (d) Spectral density 
along the high symmetry direction of the moire pattern marked with a blue line on 
panel (b). 
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on graphene [27], but it does on the clean Ru patches, forming a well known (2x2) 
superstructure epitaxially aligned with the Ru lattice, which is imaged with large 
corrugation by STM [2H]. The 0(2x2) superstructure is, thus, employed as a ruler 
to reveal the epitaxial relationship between graphene and Ru. 

Figure 3 shows an STM image of the edge of a graphene island with its 
characteristic bumps due to the moire superstructure and the simultaneously resolved 
(2x2) arrangement of oxygen atoms adsorbed on the adjacent Ru area. The green 
line follows one of the high-symmetry directions of the (2x2) superstructure, i.e. a 
symmetry direction of the Ru lattice. It deviates from the line formed by the bumps 
by 5.0° ± 0.5°. The same happens for the equivalent high symmetry directions. This 
is confirmed by the Fourier transform of the STM image shown in the inset of Figure 
3. The broad spots corresponding to the (2x2) superstructure are visible, while the 
sharper ones reflect the graphene superstructure. The moire superstructure is rotated 
with respect to the 0(2x2) pattern and, accordingly, with respect to the Ru lattice by 

^Ru,moire = 5.0° ± 0.5°. 

Due to the magnifying effect of the moire pattern [20], that is 9.96 for gr/Ru(0001) 
system, the misalignment of the C and Ru lattices can be determined with high precision 
from the observed angles between the moire superstructure and the Ru or C lattices 
and it turns out to be fgr,B.u = 0.5° ± 0.05°. This small rotation might explain the 
contradiction between the periodicity (25x25)[Il] obtained with XRD measuring the 
deformation of the last Ruthenium layers and the moire periodicity obtained from the 
STM images [5]. The superposition of the graphene and Ru(OOOl) lattices without any 
distortion and with the [1010] directions rotated by 0.5 produces a (24x24) periodicity 
when the registry with the Ru(OOOl) is taken into account. This larger periodicity 
probably reduces the compression needed to accommodate the graphene overlayer over 
Ru(OOOl) (~0.78% from DFT calculations [15]) and, presumably, also the structural 
crumpling of the graphene overlayer. 

The structural corrugation of the graphene monolayer can be determined by HAS. 
The advantages of using a beam of He atoms of thermal energy (10-100 meV) as a 
probe of the surface structure are its combination of low energy with short wavelength, 
its inert and neutral character and its large cross-section for defects. As a consequence, 
HAS is a unique non-destructive and surface sensitive technique, with high sensitivity 
to low mass atoms, such as C, or light adsorbates, such as H [29]. The actual diffraction 
grating is the periodic modulation of the repulsive part of the He-graphene potential at 
the energy of the incoming He atoms. These classical turning points define a corrugation 
function ^(x,y), which is a replica of the total surface electron density profile at about 
2-3 A above the nuclei. The amplitude of the corrugation function dictates the intensity 
of the diffraction beams. The general problem of calculating diffraction intensities 
for a given scattering geometry and corrugation function consists in solving the time- 
independent Schrodinger equation with a realistic soft potential V(r). This problem 
can be solved exactly in the most general case using the close-coupling method [30]. In 
our case, we have solved the close-coupling equations applying the procedure developed 
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Figure 3. 300x270 A^ STM image of tlie edge of a large graphene island grown 
on Ru(OOOl) coexisting with a (2x2) superstructure of oxygen chemisorbed on the Ru 
patches not covered by graphene (image setpoint Vs — -IV and It = 0.1 uA). The green 
line guides the eye along the high symmetry direction of the 0(2x2) superstructure. 
The blue line guides the eye along the moire pattern. The color scale has been adjusted 
independently on the (2x2) superstructure and on the graphene island. Inset: Fourier 
transform of the STM image. 



by Manolopoulos et al. [31], which achieves convergence much faster than the method 
originally proposed by Wolken [30] and is therefore more appropriate for calculations of 
large unit cells like the moire of graphene on Ru(OOOl). 

Figure 4(a) shows the in-plane (0=0°) (black line) and out-of-plane (0=1.3°, 15.5°) 
(green and red lines respectively) He- diffract ion spectra from a graphene overlayer grown 
on Ru(OOOl) measured along the [1120] direction, with a He beam energy of 43 meV. 
The intensity of the specularly reflected (0,0) peak depends on the structural perfection 
of graphene, reaching 25% of the incident beam for highly perfect layers. The in-plane 
He-diffraction spectrum shows both the first peaks of the moire superstructure (close 
to the specular beam) and the (1,1) of the C lattice. It is worth mentioning that the 
appearance of the first moire peaks and the (1,1) in the same in-plane scan cannot 
be taken as an indication that the corresponding lattices are aligned. The inset in 
Figure 4(a) shows the reciprocal lattices of the graphene atomic structure (black dots) 
and the moire superstructure (red dots). For the moire superstructure two reciprocal 
lattices, corresponding to two domains rotated ±5° are shown. With the existing 
angular resolution (1.4°), estimated from the full width at half maximum (FWHM) 
of the measured peaks, and the small deviation between the atomic rows and the moire 
pattern a scan through the atomic lattice peaks in the reciprocal space unavoidably 
passes through the first two peaks of the moire superstructure (note the inset in Figure 
1(b)). 
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Figure 4. (a) Angular distribution of the He atoms scattered off graphene/Ru(0001). 
The azimuthal angle, 0, is defined respect to the [1120] direction. The energy of the 
incident He atoms is 43 meV, the incident angle respect to the surface normal is 49° 
and the surface is kept at lOOK during the scattering experiments. The inset indicates 
the two rotated (±5°) moire superstructure reciprocal lattices (red hexagons) and the 
graphene atomic lattice (black hexagon), (b) Evolution of the intensity of the specular 
beam and the first diffraction peaks of the moire superstructure as a function of the 
azimuthal angle (p. As we move out-of-plane, cf) increases, and the intensity of the 
specular beam and the first order diffraction peaks decreases. For larger than 1.2°, 
the specular peak disappears, and the first out-of-plane diffraction peaks from the 
moire superstructure appear. The inset shows the in-plane scattering geometry. 



The relative intensity of the different diffraction peaks with respect to the specular 
peak determines unequivocally the corrugation function, ^(x,y), i.e. the corrugation of 
the constant charge density contour where the He atoms are reflected. The potential 
V(r) has been modeled by a discrete Fourier series, 

y (r) = Y. Ace-^^^e'^^^"'^) , (1) 
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where R is the component of r in the surface plane and G a reciprocal lattice vector. Aq, 
and i?G are the coefficients corresponding to the amplitude and exponential attenuation, 
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Figure 5. Close-coupling calculations of the in-plane diffraction probabilities for two 
different corrugation amplitudes of the moire superstructure, 0.85 A (left panel) and 
0.15 A (right panel) with an incident energy of 43 meV. Note the large differences in the 
relative values of diffraction probabilities with respect to the specular beam (0,0). The 
table in the right panel compares the best fit intensities obtained with close-coupling 
calculations with the experimental data of the moire superstructure shown in Figure 
4(a), corresponding to 6*^=49° and Ei= 43 meV. The intensities are normalized to the 
specular peak. Inset: real space potential profiles at the corner (1) and in the middle 
(2) of the moire unit cell. 
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respectively, of the different terms used in the fitting procedure (two terms in our case). 
Vg=o denotes the laterally averaged potential, and Vq^o the Fourier coefficients of the 
periodic part of the potential. Vg=o has been modelled by a two parameter Morse 
potential, 

V{z) = D(e-2°(^--o) 



2g-a{.-^o))^ 



(2) 

using for D (potential well depth) and a (range parameter) the values D = 14 meV and a 
= 1.15 A~^, which were derived from the selective adsorption resonances reported for the 
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He-graphite(OOOl) interaction by Boato et al. [32J. The unit cell was modelled including 
two different coefficients Aq, corresponding to the graphene atomic periodicity and the 
moire periodicity (but not the small misalignment between them), whereas the same 
coefficient Bq, was used for both, and the corrugation function ^(x,y) was determined 
by fitting the measured in-plane and out-of-plane diffraction intensities by means of 
a trial and error procedure. In the table show as inset in Figure 5 (right panel) we 
compare the experimental intensity of the diffraction peaks, normalized to the specular 
intensity, with our best fit. The corresponding corrugation function obtained for the 
moire superstructure periodicity has a maximum corrugation amplitude of 0.15 A, one 
order of magnitude smaller than the value theoretically predicted ^M. 115] for the aligned 
graphene overlayer. This low value was expected from the similar intensities observed 
in Figure 4(a) for the C-associated (1,1) He diffraction peak and the ones coming from 
the moire superstructure. A possible explanation for the high corrugation predicted by 
the DFT calculations is the forced alignment of the graphene lattice to the underlying 
Ru lattice imposed in the structural model used in the DFT calculations [HI [15] , not 
allowing the system to explore all the possible mechanisms to relax the atomic positions 
and producing an artificially large structural corrugation. The energy dependence of the 
corrugation amplitude can be found in the inset of Figure 5, which shows the potential 
in real space, giving a surface corrugation of 0.15 A at Ei^ = 20 meV.From this graph 
one can also estimate that the corrugation amplitude remains almost unchanged around 
0.15 A in the range of incident energies between 30-lOOmeV. 

The previous fit has only sense if the total intensity of diffracted and specular 
peaks is proportional to the scan performed with the existing resolution. This could 
be wrong if the width of the diffracted peaks would be substantially larger than the 
width of the specular peak in the azimuthal direction. In order to establish that 
the relative intensities of moire/specular peaks measured in a line scan reflect the 
true relationship between intensities, the width of the peaks has been determined by 
performing additional out-of-plane measurements. Figure 4(b) shows the evolution of 
the intensity of the specular beam and the first order diffraction peaks corresponding to 
the moire superstructure when the He detector is moved in tiny steps in the out-of-plane 
direction. When the detector is moved out-of-plane, both the specular beam and the in- 
plane diffraction peaks corresponding to the moire periodicity get smaller in such a way 
that the intensity ratio between the specular peak and the first order diffraction ones 
remains almost constant (curves measured with (j) = 0.1°, 0.4°, 0.5° and 0.6° in Figure 
4(b)). For larger (p values (1.2° and 1.6°), we observe the appearance of the first order, 
out-of-plane diffraction peaks corresponding to the moire superstructure. We explored in 
this way cp values in the range ±20°, which correspond in the reciprocal space to more 
than one unit cell of the atomic periodicity, since the first order out-of-plane peaks 
corresponding to the atomic periodicity appear at 15.5°, as can be seen in Figure 4(a) 
(red curve). These measurements allow us to conclude that the FHWM of diffraction 
peaks corresponding to the (1x1) and the moire superstructure are comparable. This 
means that there is no more intensity than the one we show in these spectra and. 
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therefore, our comparison in intensity between the specular peak and the diffraction 
ones is meaningfuL 

Finally, we have studied the sensitivity of our best-fit results to the Morse 
parameters. Briefly, we performed several fits varying D between 12-18 meV, and a 
from 0.8 to 1.4 A . As a result, we found that our data could always be fitted with a 
corrugation amplitude in the range 0.15-0.20 A. This means that with our current data, 
although we cannot refine further the best-fit parameters, we can safely state an upper 
limit to the maximum corrugation amplitude of this surface at 0.20 A. Further support 
to this conclusion comes from calculations performed for the best-fit D and a values, 
which show that the intensity of the first order moire peak strongly depends on the 
corrugation's value, reaching ~ 50% of the specular peak for a corrugation amplitude 
of just 0.20 A. In addition, for this value already the second order moire peak starts 
to be visible. Since we don't see this peak in the experiment, we can conclude that the 
corresponding corrugation amplitude cannot be larger than 0.20 A. 

We have also studied the energy dependence of the diffraction intensities by 
performing measurements at higher incident energies (up to 125 meV) which were 
complemented by close-coupling calculations. These measurements show that the 
relative intensity of the first moire peak (compared to the specular one) increases from 
ca. 12% at 30 meV to ca. 27% at 125 meV, while the intensity of the (1,1) peak remains 
almost unchanged in the same energy range. This observation already shows that the 
energy dependence of the corrugation amplitude, if any, can be only very weak. We have 
also quantified this by fitting the data at the highest energy (varying also slightly the 
parameters of the Morse potential around our best-fit values at Ej = 43 meV), obtaining 
a corrugation amplitude of 0.17 A, i.e. very similar to the value we got at lower energies. 

In order to check how sensitive are our close-coupling calculations to changes 
in the surface corrugation, we performed calculations using two different values for 
the corrugation amplitude of the moire superstructure, without changing any other 
parameter in our model. In Figure 5 (right panel) we show the result of such a calculation 
for a moire corrugation of 0.15 A. Most of the intensity (up to 44%) is concentrated 
in the specular beam, and only around 6% goes to the first in-plane diffraction peaks 
due to the moire superstructure. This calculated diffraction spectrum closely reproduce 
our in-plane measurements, as can be seen from the comparison of the normalized 
intensities in the table shown in Figure 5 (right panel). As we increase the corrugation 
of the moire superstructure more intensity is diffracted away from the specular beam 
[29] . If the corrugation of the moire superstructure is increased up to 0.82 A |T7j, the 
new calculated diffraction spectra (left panel in Figure 5) show a very weak specular 
peak with less than 0.1 % of the total diffraction intensity, and most of the intensity 
is distributed among the many diffraction peaks due to the periodicity of the moire 
superstructure. This new calculated diffraction spectra does not reproduce any of the 
main features of our experimental data, in which the most intense peak is the specular 
one and the diffraction peaks due to the moire superstructure rather weak (see Figure 
4(a)(b)). 
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Figure 6. Inversion of tlie contrast in STM images of graphene/Ru(0001). (a), (b) 
and (c) show topographic STM images at three selected sample bias voltages, (d) Local 
Tunnel Spectroscopy recorded with the tip located on top of the high (black curve) 
and low (red curve) areas of the moire pattern in panel (a) . (e) Voltage dependence 
of the apparent topographic corrugation for many different experimental conditions 
(tip, tunneling current, samples and temperature). The data clearly show a contrast 
inversion for bias voltages larger than +2.6 eV. 
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It should be mentioned that on some metalhc substrates [221 EH ES] the He 
diffraction technique underestimates the surface corrugation due to the interaction of 
the He Is electrons with the metallic surface electronic structure producing the so called 
" anticorrugation effect". Calculations suggest an upper limit for this "anticorrugation 
effect" of the order of 0.2 A [351 EE] for metallic surfaces. Without a detailed calculation 
we cannot exclude the existence of this phenomenon in this system, where the electronic 
structure is strongly modulated |6]. Therefore our estimate might be a lower bound to 
the actual corrugation and, we can safely conclude that the structural corrugation of the 
moire superstructure is in between 0.15 A and 0.4 A, still much smaller than the value 
predicted by the DFT calculations [T^ [15] or the experimental values determined by X- 
ray surface diffraction [11] or low energy electron diffraction [T8j . A possible problem in 
these diffraction techniques is the low sensitivity to light atoms (i.e. carbon) compared 
to the transition metal substrate (ruthenium) which makes difficult the precise location 
of the C atoms, specially in view of the size and complexity of the unit cell with more 
than 1000 atoms for the periodicity suggested by X-ray surface diffraction [TTj . 

The STM images measured on this surface present always the same lateral 
periodicity but depending on the bias voltage the apparent corrugation of the moire 
superstructure changes strongly with voltage and can even be inverted [261 EZ]- If the 
corrugation measured with the STM has its origin in a modulated electronic structure 
[HI [13] in addition to a physically corrugated graphene layer, the value measured with 
the STM should change dramatically with the sample bias voltage. Figure 6 (a),(b),(c) 
show STM topographic images recorded at representative sample biases and at 4.6 K. A 
defect in the superstructure encircled in blue provides with an absolute reference against 
possible drift between the images. The blue line highlights that the regions visualized 
as bumps at negative sample voltages (occupied states), are seen as depressions at +2.8 
V (empty states). Panel (e) shows the voltage dependence of the apparent topographic 
corrugation for many different experimental conditions (tip, tunneling current, samples 
and temperature). When imaging occupied states, the apparent corrugation is rather 
constant (~1 A), as expected because electrons at the Fermi level contribute the most 
to the tunneling current, but the corrugation decreases continuously when injecting 
electrons in the empty states of graphene and becomes negative above +2.6 V. This 
behavior is fully reversible and do not depend on the sample temperature in the range 
between 4.6K and 300K. 

Spatially resolved scanning tunneling spectroscopy (Figure 6 (d)) shows that the 
inversion of the contrast is due to the appearance of an intense peak in the empty states 
at +3.0 eV that is localized in the apparent depressions of the topographic images 
recorded at negative bias voltage. As demonstrated elsewhere, the origin of this state is 
the hybridization of the graphene first empty state with the Ru conduction band [38]. 
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3. Conclusions 

In conclusion, we have found for gr/Ru(0001) that the graphene lattice is rotated with 
respect to the Ru lattice, a fact overlooked in previous theoretical calculations. The 
structural corrugation of the total charge of the ripples as determined by HAS is 0.15 
A, substantially smaller than predicted by DFT calculations. Most of the apparent 
corrugation in STM images is of electronic origin, being this component so strong that 
it leads to a reversible inversion of the contrast above +2.6 V. This is originated by an 
empty electronic state derived from the Ru bands and spatially localized in the valleys 
of the structural corrugation. The electronically corrugated graphene overlayer is a 
self-nanostructured playground where new physics and spatially organized chemistry is 
bound to appear. 

4. Acknowledgments 

Financial support by the Ministerio de Educacion y Ciencia through projects 
CONSOLIDER-INGENIO 2010 on Molecular Nanoscience and FIS2007-61114 and 
Comunidad de Madrid through the program NANOBIOMAGNET S2009/MAT1726 is 
gratefully acknowledged. A. P. thanks the Ministerio de Educacion for financial support. 

References 

[1] A. K. Geim. Graphene: Status and Prospects. Science, 324:1530-1534, 2009. 

[2] A. K. Geim and K. S. Novoselov. The rise of graphene. Nat. Mat., 6:183-191, 2007. 

[3] F. Schedin, A. K. Geim, S. V. Morozov, E. W. Hill, P. Blake, M. I. Katsnelson and K. S. Novoselov. 

Detection of individual gas molecules adsorbed on graphene. Nat. Mat., 6:652-655, 2007. 
[4] Ch.-H. Park, L. Yang, Y.-W. Son, M. L. Cohen and S. G. Louie. Anisotropic behaviours of massless 

Dirac fermions in graphene under periodic potentials. Nat. Phys., 4:213-217, 2008. 
[5] S. Marchini, S. Giinther, J. Winttcrhn, Phys. Rev. B, 76:075429-1-9, 2007. 
[6] A. L. Vazquez de Parga, F. Calleja, B. Borca, M. C. G. Passeggi, J. J. Hinarejos, F. Guinea and R. 

Miranda. Periodically Rippled Graphene: Growth and Spatially resolved Electronic Structure. 

Phys. Rev. Lett., 100:056807-1-4, 2008. 
[7] P. W. Sutter, J.-I. Flege and E. A. Sutter. Epitaxial graphene on ruthenium. Nat. Mat, 7:406-411, 

2008. 
[8] A. T. N'Diaye, S. Bleikamp, P. J. Feibelman and Th. Michely. Two-Dimensional Ir cluster Lattice 

on a Graphene Moire on Ir(lll). Phys. Rev. Lett, 97:215501-1-4, 2006. 
[9] J. Coraux, A. T. N'Diaye, C. Buss and Th. Michely. Structural Coherency of Graphene on Ir(lll). 

Nano Lett., 8:565-570, 2008. 
[10] T. Fujita, W. Kobayashi and Ch. Oshima. Novel structures of carbon layers on a Pt(lll) surface. 

Surface Interface Anal, 37:120-123, 2005. 
[11] D. Martoccia, P. R. Willmott, T. Brugger, M. Bjorck, S. Giinther, C. M. Schleptz, A. Cervehino, 

S. A. Pauh, B. D. Patterson, S. Marchini, J. Wintterlin, W. Moritz and T. Greber. Graphene 

on Ru(OOOl): A 25x25 Superceh. Phys. Rev. Lett, 101:126102-1-4, 2008. 
[12] A. B. Preobajenski, M. L. Ng, A. S. Vinogradov and N. Matersson. Controlling graphene 

corrugation on lattice-mismatched substrates. Phys. Rev. B, 78:073401-1-4, 2008. 
[13] A. L. Vazquez de Parga, F. Calleja, B. Borca, M. C. G. Passeggi, J. J. Hinarejos, F. Guinea and 

R. Miranda. Vazquez de Parga et al. Reply, Phys. Rev. Lett, 101:099704-1-1, 2008. 



Electronic and Geometric Corrugation of Periodically Rippled, Self-nanostructured Graphene Epitaxially Gr 

[14] B. Wang, M.-L. Bocquet, S. Marchini, S. Giinther and J. Wintterlin. Chemical origin of a graphene 

moire overlayer on Ru(OOOl). J. Phys. Chem. Chem. Phys., 10:3530-3534, 2008. 
[15] D. Jang, M.-H. Du, S. Dai. First principles study of the graphene/Ru(0001) interface. J. Chem. 

Phys., 130:074705-1-5, 2009. 
[16] X. Peng and R. Ahuja. Epitaxial graphene monolayer and bilayers on Ru(OOOl): Ab initio 

calculations. Phys. Rev. B, 82:045425-1-5. 
[17] D. Martoccia, M. Bjorck, CM. Schleptz, T. Brugger, S.A. Pauh, B.D. Patterson, T. Greber 

and P. R. Willmott. Graphene on Ru(OOOl): A corrugated and chiral structure. New J. Phys, 

12:043028-1-12, 2010. 
[18] W. Moriz, B. Wang, M.-L. Bocquet, T. Brugger, T. Greber, J. Wintterlin and S. Giinther. Structure 

determination of the coincidence phase of graphene on Ru(OOOl). Phys. Rev. Lett., 104:136102- 

1-4, 2010. 
[19] P. Sutter, M. S. Hybertsen, J. T. Sadowski and E. Sutter. Electronic Structure of Few-Layer 

Epitaxial Graphene on Ru(OOOl). Nano Lett., 9:2654-2660, 2009. 
[20] A. T. N'Diaye, J. Coraux, T. N. Plasa, C. Bussc and Th. Michely. Structure of epitaxial graphene 

on Ir(lll). New J. Phys., 10:043033-1-15, 2008. 
[21] M. Gsell, P. Jakob and D. Mcnzcl. Effect of Substrate Strain on Adsoprtion. Science, 280:717-720, 

1998. 
[22] F. Calleja, V. M. Garcia-Suarez, J. J. Hinarejos, J. Ferrer, A. L. Vazquez de Parga and R. Miranda. 

Relationship between strain and the surface electronic structure of Cu(lll) films on Ru(OOOl): 

Theory and experiment. Phys. Rev. B, 71:125412-1-12, 2005. 
[23] L Horcas, R. Fernandez, J.M. Gomez- Rodriguez, J. Colchero, J. Gomez-Hcrrero and A. M. Baro. 

WSXM: A software for scanning probe microscopy and a tool for nanotcchnology. Rev. Sci. 

Instrum. 78:013705-1-8, 2007. 
[24] E. Loginova, N. C. Bartelt, P. J. Feibelman and M. F. McCarty. Factors influencing graphene 

growth on metal surfaces. New J. Phys., 11:063046-1-18, 2009. 
[25] M. Gao, Y. Pan, Ch. Zhang, H. Hu, R. Yang, H. Lu, J. Cai, S. Du, F. Liu and H.-J. Gao. Tunable 

properties of epitaxial graphene on metal substrates. Appl. Phys. Lett. 96:053109-1-3, 2010. 
[26] Y. Pan, H. Zhang, D. Shi, J. Sun, S. Du, F. Liu, H.-J. Gao. Highly ordered, millimeter-scale, 

continous single-crystalline graphene monolayer formed on Ru(OOOl). Adv. Mater., 21:2777-2780, 

2009. 
[27] B. Borca, F. Calleja, J. J. Hinarejos, A. L. Vazquez dc Parga, R Miranda. Reactivity of periodically 

rippled graphene grown on Ru(OOOl). J. Phys. Condens. Matter, 21:134002-1-6, 2009. 
[28] F. Calleja, A. Arnau, J. J. Hinarejos, A.L. Vazquez de Parga, W. Hofer, P. M. Echenique and R. 

Miranda. Contrast Reversal and Shape Changes of Atomic Adsorbates Measured with Scanning 

Tunneling Microscopy. Phys. Rev. Lett, 92:206101-1-4, 2004. 
[29] D. Farias and K.H. Rieder. Atomic beam diffraction from solid surfaces. Rep. Prog. Phys., 61:1575- 

1664, 1998. 
[30] G. Wolken Jr.. Theoretical studies of atom-solid elastic scattering: He-|-LiF. J. Chem. Phys., 

58:3047-3064, 1973. 
[31] D. E. Manolopoulos, R.E. Wyatt and D.C. Clary. Iterative solution in quantum scattering theory. 

The log derivative Kohn approach. J. Chem. Sac. Faraday Trans., 86:1641-1648, 1990. 
[32] G. Boato, O. Cantini and R. Taratek. Study of Gas-Graphite Potential by Means of Helium Atom 

Diffraction. Phys. Rev. Lett., 40:887-889, 1978. 
[33] K. H. Rieder and N. Garcia. Energy Dependence and Softness of the Potential for He Scattering 

from Ni(llO). Phys. Rev. Lett, 49:43-46, 1982. 
[34] K. H. Rieder, G. Parschau and B. Burg. Experimental evidence for anticorrugation effects in 

He-metal interaction at surfaces. Phys. Rev. Lett., 71:1059-1062, 1993. 
[35] M. Petersen, S. Wilke, P. Ruggerone, B. Kohler and M. Schcfflcr. Scattering of Rare-Gas Atoms 

at a Metal Surface: Evidence of Anticorrugation of the Helium- Atom Potential Energy Surface 

and the Surface Electron Density. Phys. Rev. Lett, 76: 995-998, 1996. 



Electronic and Geometric Corrugation of Periodically Rippled, Self-nanostructured Graphene Epitaxially Gr 

[36] J.F. Annett and R. Haydock. Hybridization interaction between helium and a metal surface. Phys. 

Rev. B, 34:6860-6868, 1986. 
[37] B. Borca, S. Barja, M. Garnica, J.J. Hinarejos, A.L. Vazquez de Parga, R. Miranda and F. Guinea. 

Periodically modulated geometric and electronic structure of graphene on Ru(OOOl). Semicond. 

Set. TechnoL, 25:034001-1-7, 2010. 
[38] B. Borca, S. Barja, M. Garnica, D. Sanchez-Portal, V.M. Silkin, E.V. Chulkov, F. Hermanns, J. 

J. Hinarejos, A. L. Vazquez de Parga, A. Arnau, P. M. Echenique and R. Miranda. Potential 

Energy Landscape for Hot Electrons in Periodically Nanostructured Graphene. Phys. Rev. Lett., 

105:036804-1-4, 2010. 



